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1.0 Introduction 

All the equations used by these codes, expect for an unequal pipe 
split, are developed in UAH-RI reports 5-32176 and 5-32441. These 
rotes will not present all the equations used. There is a Users 
Manual available for each of these codes. 

The piping in a liquid rocket can assume complex configurations due to 
multiple tanks, multiple engines, and structures that must be piped 
around. The capability to handle some of these complex configurations 

have been incorporated into the Low Frequency Code (ADMIT) , the 

Nyquist Code (NYQUIST) , and the Intermediate Frequency Code (SSFREQ) . 
The capability to modify the input on line has been implemented. 

The configurations allowed include multiple tanks, multiple engines, 
the splitting of a pipe into unequal segments going to different (or 
the same) engines. Straight pipes, bends, inline accumulators, tuned 
stubs, Helmholtz resonators, parallel resonators, and pumps are the 
types of piping that may be used. 

The three codes require the computation of admittance looking toward 
the tank. Therefore, they all have a number of similar routines for 
reading in the tank, piping, and engine descriptions. They also have 
similar routines for computing the admittance looking toward the tank. 
ADMIT uses only one type of feedline, either a LOX or FUEL line. 
NYQUIST can use no feedlines, one feedline (either LOX or FUEL), or 
two feedlines (both LOX and FUEL) . SSFREQ requires both feedlines to 
be present (LOX and FUEL) . 

The common features of the three codes are discussed first, with two 
sample runs described. Then each of the three codes are discussed 
along with the sample runs. Several on-line variations are suggested 
for the second sample run. Input descriptions and the input for the 
sample runs are given in the appendices. 

Homework assignments for the codes are given in appendix 
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2.o Common Features of the Codes 


All three programs compute the admittance looking toward the tank and 
uses this admittance for other calculations. The tank, pa.pi ng , 
engine parameters are used in this calculation. The following types 
of piping may be used: straight sections, bends, inline accumulator or 
manifold, tuned stub, Helmholtz resonator, parallel resonator, pump. 

In addition, the pipe may be split into several pipes, either 
identical or different. If there is a split in the PiPJ?*' a “ 
iteration is required to compute these impedances. Multiple tanks and 

engines may be used. 


2.1 Sample Runs 

The codes should be run from hard disk because the use temporary files 
which may exceed the free space on the floppy disk containing the 
codes. The codes are started by typing the name of the code (ADMIT, 
NYOUIST, or SSFREQ) . After an introductory message, you are asked 
what units to use for frequency (rad/sec or Hertz). Then you * re 
asked questions about the input files. If a split is present in the 
piping, a default value of 20 iterations is set. Then you are asked 
if you wish to change the number of iterations. If you are making 
some quick and dirty runs, lower the number; if making a final run, 
you may want to increase the number. 


After the engine and piping files are read, questions about keyboard 
input or special files are asked. Where more than one value is 
requested, the input is in free format. That is, the values may be 
separated by spaces, commas, or carriage returns. For example, the 
following inputs give identical results: 

Case #1. 1 40 40 

Case #2. 1,40,40 

Case #3. 1 

40 

40. 

Note that real values may be given with or without a decimal point 
if it is an integral value. Be sure to include 0's in file input. 


After the computations are complete, you are prompted for the 
graphical output desired. After a graph is made, the program waits 
for you to hit a carriage return <CR> before continuing. After the 
plots are finished, you are asked what you would like to do next, 
end, run a new case, or possibly run a new frequency range with the 
old data. At various points in the run you will be given a chance 
to modify data, rewind an input file, or bypass the current data. 
When a frequency range is asked for, if 0 0 0 is input the code will 
bypass that set of data. 

Remember, graphics must have been run before running the application 
code to be able to dump the plots to the printer. 


TAere are two sample runs included on the disks. The first is a 
complex configuration with two LOX and two FUEL tanks feeding three 
engines. One line splits into two unequal pipes going to two 


- 2 - 



different engines. The second tank feeds one engine. The second 
sample run is a straight pipe which will be used for the on-line 
variations . 


The first sample run is to give the student a chance to run a complex 
problem. This problem will introduce the student to the code, keyboard 
input required, the graphs produced, and the time required to run a 
complex problem. The following table shows the run time (after 
keyboard input) on several computers for this problem. The 286 and 
386 had co-processors. 


Program 

286/8 

time / factor 

386SX/16 
time / factor 

486/33 

time / factor 

ADMIT 

133.03 / 26.14 

32.41 / 6.37 

5.09 / 1.00 

NYQUIST 

684.43 / 25.71 

165.55 / 6.22 

26.62 / 1.00 

SSFREQ 

1041.06 / 28.08 

242.72 / 6.54 

37.08 / 1.00 


Where time is in seconds and factor is time divided by the 486/33 
time. The runs made were for the following: 

Hertz ENG.RLN & LOX.RLN 
frequency range: 1,40,40 

pipe sections: 5,9,10,4,5,7,7,7,4,3,3,4,3,3,4,3,4,3,4, 

3,4,5,9,10,4,5,7,7,7,4,3,3,4,3,3,4,3,4 

Hertz ENG.RLN, LOX.RLN, & CONST. RLN (no fuel) 
frequency range: 1,40,300 

Hertz ENG.RLN, LOX.RLN, FUEL. RLN, & IMODE.RLN 
frequency range: 110,116,3 
tau range: . 0005, . 0025, 11 

The second sample run is a simple configuration which is used to guide 
the student in on-line modification of data and the effects of certain 
variables on the output. Changes in geometry will be made in ADMIT 
and NYQUIST. The input unique to NYQUIST and SSFREQ will also be 
varied. When making variations, rewind those files not involved in 
the changes. This should help the student become familiar with the 
codes and also the effect of certain input on results. 


ADMIT 

NYQUIST 

SSFREQ 


3.0 Program ADMIT 

Program ADMIT uses the impedance (reciprocal of admittance) looking 
toward the tank and also the impedance looking toward the engine in 
the calculation of the pressure transfer function. This is the only 
one of the programs which will plot the admittance looking toward the 
tank. The piping layout is plotted on the same graph as the 
admittance. Two other plots are available: a 3— D surface plot of the 
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pressure transfer function and a contour plot of the 
function. The plots may be in black and white or in 
user's choice. 


pressure transfer 
the color of the 


3.1 Sample Runs 

ADMIT displays the maximum number of frequencies and maximum number 
of (integration) points along the pipe. These values are a function 
of the amount of memory available for programs. If these 
too small for your problem, try getting rid of most of the TSR s 
loaded. This will probably require changing the AUTOEXEC.BAT and 
possibly the CONFIG.SYS files and restarting the computer. 

The maximum number of points along a pipe is measured from a tank to 
an engine. For sample run # l, there are three pipes: tank #1 to 
engine #1, tank #1 to engine #2, and tank #2 to engine #3. Each of 
these pipes may contain the maximum number of points. The number of 
points is the cumulative total of points en each section from the tank 
to engine. Thus, a split pipe has the same number of points from the 
tank to the split, but, may have a different number of points from the 
split to the engine. 

The admittance calculation is not affected by the number of points 
along a pipe, only the pressure transfer function. Therefore, if you 
are only interested in the admittance, set the number of points for 
each pipe section to one. If you are interested in the pressure 
transfer function, use near the maximum. You may try to set the 
number of points of a section to the length of the section. If this 
gives too many points, try half that value. Keep reducing the numb 
of points until the total is less than the maximum. Try to keep at 
least two points per pipe section. 


Sample run #1 using files ENG.RLN and LOX.RLN 


Run a frequency range of 1 40 40. Make admittance plots of 
lines. Change colors of the plot. Then make surface plots 
various options. Likewise, do contour plots of the lines, 
session. 


the three 
trying the 
End the 


Sample run #2 using files CLASS. ENG and CLASS. LOX 

Start program by typing ADMIT and specifying Hertz. Remember to make 
the input files CLASS. ENG and CLASS. LOX. Again make the frequency 
range 1 40 40. This time, simply make the number of pipe segments 
equal to the length of the pipe (e.g. the first pipe is 20 feet long, 
therefore break it into 20 segments). Now look at the results. 


3.2 On-line Variations 

After looking at the results, change the piping for variation 1 as 
follows: 


section 1 
section 2 


leave as is 
leave as is 
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section 3 
section 4 
section 5 
section 6 
section 7 
section 8 
section 9 


leave as is 
leave as is 

change to st. section 50 

leave as is 

leave as is 

leave as is 

leave as is 


ft. long & 1 ft. diameter 


After looking at the results, change the piping for variation 2 as 
follows: 


section 1 
section 2 
section 3 
section 4 
section 5 
section 6 
section 7 
section 8 
section 9 


leave as is 
leave as is 
leave as is 

leave as is _ 

change to st. section 30 ft. long & 2 ft. diameter 

leave as is 
leave as is 
leave as is 
leave as is 


After looking at the results, change the piping for variation 3 as 
follows: 


section l 
section 2 

section 3 
section 4 

section 5 
section 6 

section 7 
section 8 

section 9 


leave as is 

change to a -90 deg. bend with 1.16245 ft. radius, 
1 ft. diameter, and 0 ft. length beyond bend 
leave as is 

change to a 90 deg. bend with 1.16245 ft. radius, 

1 ft. diameter, and 0 ft. length beyond bend 
change to st. section 30 ft. long & 1 ft. diameter 
change to a 90 deg. bend with 1.16245 ft. radius, 

1 ft. diameter, and 0 ft. length beyond bend 
leave as is 

change to a —90 deg* bend with 1*16245 ft* radius, 
1 ft. diameter, and 0 ft. length beyond bend 
leave as is 


Now look at the results. Note that the admittance curve changed with 
the first two variations, but changed very little with the addition of 
the bends (that kept the length and diameter of the piping the same as 
the straight pipe) . The pressure transfer function changed 
considerably for all the three variations. If the admittance is 
relatively unchanged, the Nyquist and n-T plots are unchanged. 


4.0 Program NYQUIST 

The Nyquist program can use the admittances Gox and Gf. The code was 
written to plot the Nyquist curves for the four cases: neither 
admittance used, Gox only, Gf only, and both admittances used. Three 
types of plots are available: the piping layout if either or both of 
the lines are present, Nyquist curve in the complex plane, and the 
phase-gain plots of the Nyquist equation. 
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in addition to the Nyquist plots of these four equations, Phase-Gain 
plots are also available. 


4.1 Sample Puns 

S am ple run #1 using files ENG.RLN, LOX.RLN, and CONST. RLN 

Run with only the LOX data with a frequency range of 1 40 300. Plot 
the three LOX lines* After computations are complete, try the plot 
available. When through looking at the curves, exit the program. 


Sample run #2 using files CLASS. ENG, CLASS. LOX, and CLASS. CON 

Again run only the LOX data with a frequency range of 1 40 300. When 
finished, make the following variations. 


4.2 On-line Variations 

For variations, make the same changes as for ADMIT (change length, 
diameter, and add bends. Then, change the transport lag (TAUT) by 
a factor of two. 


5.0 Program SSFREQ 

Proqram SSFREQ works in the intermediate frequency range where the 

piping, engine, and nozzle all interact. The program computes n 

(pressure interaction index) for a range of frequencies and a range of 
r * s (sensitive time lag). After the range of t's for a given 
frequency have been run and the n's displayed on the screen, the user 
may request a plot of n vs t for that frequency. After the range of 
frequencies have been run, n vs r is plotted for all frequencies on 
one graph. 

5.1 Sample Runs 

Sample run #1 using files ENG.RLN, FUEL. RLN, LOX.RLN, and IMODE.RLN 

For this run, you may enter the frequency range (110 116 3) and the 
tau range (.0005 .0025 11) or read them in with files IMODE.FRQ and 
IMODE.TAU. After each set of taus for a frequency is computed, you 
may wish to see the n-T curve for that frequency. Simply enter y to 
see the curve or no to bypass the plot. After the complete frequency 
range is computed, a composite plot with all frequencies is presented 

for each engine. 

Sample run #2 using files CLASS. ENG, CLASS. FUL, CLASS. LOX, and 

CLASS . IMO 

Start as for sample run #1 except change all the input files to those 
starting with CLASS. After looking at the results, begin the changes. 
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5.2 On-line Variations 

SSFREQ does not allow changes to the piping. However, the data 
unique to SSFREQ may be changed. Try the following variation . 

1 . increase diameter (CDIAM) by 20% 

2. increase length (XLCD) by 10% 

3. increase liquid velocity (ULOD) by 20% 


6.0 Equations 


6.1 Piping Equations 

In the following equations, n = s/a. 

1. Bend - The effective straight pipe length and diam ®*® r a *® computed 
using the equations on page II-4 of report UA J *“*?*’? ht 
(curves on page II-II have been curve fit) . The g 

pipe equations are then applied to the effective length and 

diameter. 


2. Straight Pipe 

Zt(I-D + Z 0 (I) -tanh(n-l) ^ 

Z t (I) = Z 0 (I) * C — 1 

* Zq (I) + Z t (I-l)-tanh(n.l) 

Z g (I) = {e n * *1 • [Z 0 (I+l) + z g (I+l)]-(l - N-M-e” 2,11,1 ) 

- Z 0 (I+1) • (1 - N* e” 2 * n * l*e 2 * n * ^1) >/ (1 + N* e" 2 *n* 1. e 2 * n * ^1) 

where N = [Zg(I+l) - Z-t(I-l) ]/[Zo(I+l) + Z^(I-l)] 

M = [Z 0 (I+l) “ Zg(I+l) ]/[Z 0 (I+l) + Zg(I+l)] 

1 = L ( I ) + L (I+l) 

lj. = L (I+l) 

3. Inline Accumulator 
Z e =l/(C-s) 

Z t (I) = Z e *Z t (I“D / [Zt + 2 e ] 

Z g (I) = Z e -Zg(I + l) /[Zg(I + D + Z e ] 

4. Tuned Stub 

Z e = Zo/tanh(n*l) 

Z t (I) = Z e -Z t (I-l)/[Z t (I-l) + z e ] 

Zg ( I ) = Z e *Zg(I+l)/tZg(I+l) + Z e ] 


- 7 - 



5. Helmholtz Resonator 

Z e = (1 + L*C* s 2 ) / (C* s) 

Zt(I) = Z e -Z t (I-l)/[Zt( I - 1 > +z e] 

Zg(I) = Z e -Zg(I+l)/[Zg(I+l)+Z e ] 

6. Parallel Resonator 
Z e = L-s/(l + L* C • s 2 ) 

Z t (I) = Zt(I-l) + z e 

Z g (I) = Zg(I+l) + Z e 

7. Pump 


3 P 



Z t (I) = {Zt (I“D + (Zp + !*• s) • [1 + Z t (I-l).Cs]>/[l + Z t (I-l).Cs] 
Zg(I) = [L • S - Zp + Zg(l+1)]/{1 + C*S*[L*S - Zp + Zg(I + l)]> 

8. Split Pipe 

The equations will be shown for a split into 3 unequal pipe segments. 
Case I. Finding the impedance looking toward the tank (Zt) • 


1-1 


z t (i-D 

Zg(I-l) 


| II | 11+1 


12 | 12+1 


13 I 13+1 


Z t (1-1) * Z t (I) 
Zg(I-l)* Zg ( I ) 


Section II, looking toward the tank sees Zt(I-l) and (m-1) Zg(I-l)*'s 
in parallel. Therefore the effective Z e (Il) is 


1 1 1 1 

= + + 

Z e (ID Z g (I2-l)* Zg(I3-l)* Zt(I-l) 

This Z e (Il) is used in the equations developed for Zt(H) instead of 
Z t (I-l) • 
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Case II. 


Finding the impedance looking toward the engine (2 g ) . 



1 11+1 

| 11+2 

1 

| 12+1 

| 12+2 

1 


13+1 


13+2 


Z t (D Z t (I)* *t< 1+1 > 

Zg(I) Zg(I) Zg(I+l) 

Section I, looking toward the engine sees m sections 1+1 in parallel, 
Therefore the effective z g (I) is a function of three pipes. Thus, 
compute Zg* for each pip® and then 


z g (l) ~ Zg(Il)* Z g ( 12 ) * Z g (13) * 

6.2 Pressure Transfer Function 

The equation for the pressure at any point in a pipe is derived on 
page 25 of NASA Contractor Report 5-32176. 


p(x,s) 

Pg(s) 

where 

where 


= ( 


Z Q e“ n ' x - N*e” n * (2‘i-x) 

) • [ r 3 


Zq + Z g 
n = s/a 

Zq “ z t 


1 - N*M* e“ 2 * n * 1 


N = 


M = 


Z 0 + z t 

Zq “ Z g 
Z 0 + Z g 


6.3 Nyquist Equations 

On page 47 of the report the following equation is derived 

e -T«s (l + z) dc* r* (1 + r) dc* 

{[1 + ( )]*Gox + I 1 : < ) ] ' G f } 


= -l, 


( l+©c * ® ) 


dr 


dr 
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in order to simplify the notation, the following definitions are used 


e -T-s 

Ki = ; 

(l+0 c -s) 


Ai 


[1 + 


(1 + r) 
C* 


3C* 

( )] 

dr 


Thus, 

The equations used are 

K ( jco) = 2-Ki 

K(jw,G ox ) = Kx*Ai 

K(jw,Gf) = Ki*A2 

K(j«,G ox ,Gf) = Kx* (A 1 +A 2 ) 

In addition to the Nyquist plots 
plots are also available. 


= - 1 . 

neither admittance used, 

G ox used, 

Gf used, 

both admittances used, 
of these four equations, Phase-Gain 


r* (1 + r) 3c* 

A 2 = [1 • <— )] 

C dr 


the equation may be expressed as Kx*^ A l* G ox + A 2’ G f) 


6.4 Intermediate Frequency Equations (n - t) 

The equations to analyze intermediate mode instabilities were 
developed in UAH report 5-32176. In some cases with long and 
laborious mathematics. The pertinent equations are presented as used 
by the code SSFREQ. 


Xi = (y - 1 ) uuq + (1 + r ) (dh^/dr) (M/s) e" ST T (G ox - rGf)Poo 
Yx = -up 0 

Zi = (l/7>up 0 + PL 8 U L 0 
Wx = 2 UUq 

Mx = M{e“ ST T[K(s) + H(S)]P 00 " P(S)P0> 

PO = POO cosh (sx) 

U 0 = “ { 1/7) Pqo Sinh(sx) 
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rx 


{[s(Wi - Xi) + M 1 ]cosh[s(x - X')] 
0 


Ui = Yi + 

+ s(Yi + Z x ) sinh[s (x - x')]}dx' 

*x 

Pl/7 = ~ W 1 " <[s(Wi - X x ) + MjJsinhtsU - x')] 

J 0 

+ s (Y x + Z^) cosh[s (x - x')]}dx' 
u' (L c ) + A p' (L c ) + C a' (L c ) = 0 


o' = a 0 + a x = -(l/u) 


{((7 - l)/7>Po + (1 + r) (dh^/dr) (G ox 
0 


- rG f )P 0 oe” ST T}M exp [-s 


rx 

(l/u)dx' ' ]dx' 
x' 


where u' = u<> + u^ and p' = Po + Pi» 

H(s) = (1 + r) [ (r/c*) (3c*/3r) - n r sf] (G ox - rG f )/r 
K (S) = (1 + ST T ) (G ox + Gf) 

P(s) = n (1 - e sT ) 

For small perturbations/ A and C may be approximated by 

A * M(y - 1) / (2-y) where M is Mach number at intersection of 

the chamber and nozzle 

C » 0 
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6.5 Nomenclature 


a - speed of sound ft/sec 

A - nozzle pressure admittance coefficient 

C - capacitance sec 

C* - characteristic velocity _ ft/sec 

C - nozzle entropy admittance coefficient 

G - admittance 

h - enthalpy (ft/sec) *2 

H - combustion response function for mixture ratio oscillations 

K - combustion response function for mass flow oscillations 


1 - length 

L - length 

L - inductance 

m - mass flow 

H - combustion response 

M - Mach number 

n - pressure interaction index 

n r - mixture ratio interaction index 

p - pressure 

Pq 0 - maximum overpressure 

P - combustion response function for pressure 
r - mixture ratio 
s - complex frequency 
u - velocity 
x - distance along pipe 
Z - impedance 

y - ratio of specific heats 
@ - characteristic time constant 
p - density 
a - entropy 
t - time lag 


ft 

ft 

sec 

lbm/sec 


lbf/ft‘2 

lbf/ft“2 


1/sec 

ft/sec 

ft 


sec 

lbm/ft~3 

sec 


Subscripts 


c - chamber 
e - effective 
f - fuel 

g - looking toward engine 
g - generator 
L - liquid 
ox - oxidizer 
p - pump 

t - looking toward tank 
T - total 
0 - lossless line 

0 - first element of expansion 

1 - second term of expansion 


Superscripts 

* - seen by split 

• - perturbation value 

° - value at ejector face 
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7.0 summary 


Three codes for the PC have been written using the admittance looking 
from the engine toward the tank. Because of the complexity of the 
codes, they should be run on a 386 with co-processor or more powerful 
PC. They will all work on most any graphics monitor (monochrome, CGA, 
EGA, or VGA) . The colors of most of the graphs may be set by the 
user. If the DOS graphics TSR has been run (i.e. typing GRAPHICS at 
the DOS prompt) , the graphs may be sent to a printer by hitting the 
PrtSc key. 

ADMIT computes the pressure tra n sfer function. The admittance curve, 
pressure transfer function surface plot (3-D) , and contour plot may 
made. Colors may be specified for each of the plots separately. The 
surface plot allows the user to change view angle, type of plot (wire 
frame or solid) , and color. Input data may be modified at the 
keyboard . 

NYQUIST solves the Nyquist equation for no lines, FUEL line only, LOX 
line only, and/or FUEL and LOX lines. Nyquist plots or phase-gain 
plots may be made for any condition run. Also, a plot of the lines 
available may be made on request. 

SSFREQ solves the intermediate frequency equations to generate a n-T 
curve. This code couples the lines, chamber, and nozzle to obtain 
the effect in the intermediate frequency range. 

All of the codes were written in Microsoft FORTRAN 5.0. ADMIT and 
SSFREQ used GRAFMATIC plot routines by Microcompatibles, Inc. The 
plot routines in FORTRAN 5.0 was used in NYQUIST. 



Appendix I 


1.1 Common Input 
Description of file ENG.RLN 



Card # 1 

MENG - number of engines 

Card # 2 

TFLOW(I) - total flow in engine 
PCHAMB(I) - chamber pressure 
DPROR(I) - pressure drop across orifice 
Read card # 2 "number of engines" times 

Description of files LOX.RLN or FUEL.RLN 

Card # 1 

TITL - LOX or FUEL title 

Card # 2 

MTANK - number of tanks 

Card # 3 

VOL (I) - volume 

LFLOW(I) - mass flow 

KTANK(I) - bulk modulus 

DENS (I) - density 

Read card # 3 "number of tanks" times 
Card # 4 

MLINE - number of lines leaving tank 

Card # 5 

ITANK(I) - tank number 

IENG(I) - engine number (0 if split follows) 
Card # 6 

SEGHN(I) - number of segments 
SPLIT (I ) - n umb er of unique splits 

Card # 7 

SECTN(I,M) - section type 
PIPE1 (I/M) - radius 
length 
or volume 

PIPE2 (I/M) - angle 

diameter 
or bulk modulus 
PIPE3 (I/M) - diameter 
volume 
dp /dm 
or 0.0 

PIPE4 ( I , M) - end length 
L 

or 0.0 

PIPE5 (I/M) - C 

or o . o 

Read card # 7 "number of segments" times 


lbm/sec 

lbf/ft~2 

lbf/ft'2 


ft~3 

lbm/sec 

lbf/ft‘2 

lbm/ft‘3 


ft 

ft 

ft~3 

deg 

ft 

lbf/ft“2 

ft 

ft* 3 


ft 

sec 


sec 
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if split > o 
Card # 8 

SEGMN(I) - number of segments 

NOLINE (I) - number of identical lines 

IENG(I) - engine number 

Card # 9 


SECTN (I/M) - section type 

PIPE1 (I/M) - radius ft 

length ft^ 

or volume ft 3 

PIPE2 (I/M) - angle de, 3 

diameter ft 

or bulk modulus lbf/ft“2 

PIPE3 (I,M) - diameter ft^ 

volume ft“3 

dp /dm 
or 0.0 

PIPE4 (I/M) - end length ft 

L sec 

or 0.0 

PIPES (I/M) - C sec 

or 0.0 


Read card # 9 "number of segments" times 
Read card # 8-9 "number of splits" times 
Read card # 5-9 "number of lines" times 


type 

name 

PIPE1 

PIPE2 

PIPE3 

PIPE4 

PIPE5 

0 

bend 

radius 

angle 

diameter 

end len. 

0.0 

1 

straight 

length 

diameter 

0.0 

0.0 

0.0 

2 

inline 

length 

diameter 

0.0 

0.0 

0.0 

3 

tuned 

length 

diameter 

0.0 

0.0 

0.0 

4 

Helmholtz 

length 

diameter 

volume 

0.0 


5 

parallel 

length 

diameter 

volume 

0.0 

warn 

6 

pump 

length 

diameter 

dp /dm 

L 

c 

7 

manifold 

volume 

bulk mod. 

o 

• 

o 

0.0 

O 

• 

o 


Dimensions: 

radius/ length/ diameter/ end length 
angle 

volume o 

dp/dm (non-dimensionalized by m/p c ) 
L 
C 

bulk modulus 


ft 

deg 

ft~3 

sec 

sec 

lbf/ft~2 
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1.2 Input Unique to NYQUIST 
Description of file CONST. RLN 


Card # 1 

TAUT (I) - transport lag 0 

CSTAR ( I ) - characteristic velocity (P c * A t/ m ) 

RBAR(I) - mixture ratio 0 

THETAC(I) - characteristic time constant (p C ‘ v c/ m ) 
DCDR(I) - dc*/dr 

Read card # 1 "number of engines" times 


sec 

ft/sec 

sec 

ft/sec 



1.3 Input Unique to SSFREQ 
Description of file IMODE.RLN 


Card # 1 

TITLF - title of IMODE file 

Card # 2 

NVAL - number of x stations 

Card # 3 

XBARD(I) - x location 

PBAR(I) - pressure 

TBAR(I) - temperature 

Read card # 3 "number of x stations" times 

Card # 4 

DTAUD - invariant time lag 

NRD - mixture ratio interaction index 

LAMDAD - damping part of frequency 

Card # 5 

CDIAM - chamber diameter 

TDIAM - throat diameter 

XLCD - length of combustion chamber 

Card # 6 

6AMMAD - ratio of specific heats 

rgas - gas constant 

POOD - maximum overpressure 

RBADR - mixture ratio 


ft 

lbf/ft~2 

°R 

sec 


ft 

ft 

ft 


(ft/sec) “2/ °R 
lbf/ft‘2 


Card # 7 

DCSDRD - dc*/dr ft/sec 

DHLDRD - dh/dr (ft/sec) ~2 

RHOLOD - mass of liquid per unit chamber volume lbm/ft~3 

ULOD - axial component of liquid velocity ft/sec 

Read card # 2-7 "number of engines" times 
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Appendix II 


II. 1 Files for First Sample Run 


Common Input for First Sample Run 
ENG.RLN file for First Sample Run: 


3 

853.5 

1707.0 

3414.0 


4 . 502040E+05 
4 . 502 04 0E+05 
4 . 502040E+05 


1.610532E+06 
1. 610532 E+ 06 
1.610532E+06 


LOX.RLN file for First Sample Run: 


Sample Run 
2 

1.956300E+04 2928.0 1.185883E+07 71.4 

1.956300E+04 2928.0 1.185883E+07 71.4 


2 

1 0 
13 2 


1 

15.0 

1.416 

0.0 

0.0 

0.0 

0 

35.0 

45.0 

1.416 

0.0 

0.0 

1 

30.0 

1.416 

0.0 

0.0 

0.0 

0 

3.5 

135.0 

1.416 

0.0 

0.0 

1 

15.0 

1.416 

0.0 

0.0 

0.0 

1 

20.641 

1.416 

0.0 

0.0 

0.0 

1 

20.558 

1.416 

0.0 

0.0 

0.0 

1 

20.558 

1.416 

0.0 

0.0 

0.0 

1 

8.541 

1.416 

0.0 

0.0 

0.0 

1 

6.383 

1.416 

0.0 

0.0 

0.0 

0 

4.25 

90.0 

1.416 

0.0 

0.0 

1 

9.33 

1.416 

0.0 

0.0 

0.0 

0 

3.33 

80.0 

1.416 

0.0 

0.0 

5 

1 1 





1 

3.53 

0.708 

0.0 

0.0 

0.0 

1 

12.2 

0.708 

0.0 

0.0 

0.0 

0 

1.28 

35.0 

0.708 

0.0 

0.0 

1 

12.2 

0.708 

0.0 

0.0 

0.0 

7 

13.5 1. 

183346E+07 0.0 

0.0 

0.0 

5 

1 2 





1 

3.53 

1.00126 

0.0 

0.0 

0.0 

1 

12.2 

1.00126 

0.0 

0.0 

0.0 

0 

1.28 

35.0 

1.00126 

0.0 

0.0 

1 

12.2 

1.00126 

0.0 

0.0 

0.0 

7 

13.5 1. 

183346E+07 0.0 

0.0 

0.0 

2 3 






18 

0 





1 

15.0 

1.416 

0.0 

0.0 

0.0 

0 

35.0 

45.0 

1.416 

0.0 

0.0 

1 

30.0 

1.416 

0.0 

0.0 

0.0 
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0 

3.5 

135.0 

1.416 

1 

15.0 

1.416 

0.0 

1 

20.641 

1.416 

0.0 

1 

20.558 

1.416 

0.0 

1 

20.558 

1.416 

0.0 

1 

8.541 

1.416 

0.0 

1 

6.383 

1.416 

0.0 

0 

4.25 

90.0 

1.416 

1 

9.33 

1.416 

0.0 

0 

3.33 

80.0 

1.416 

1 

3.53 

1.416 

0.0 

1 

12.2 

1.416 

0.0 

0 

1.28 

35.0 

1.416 

1 

12.2 

1.416 

0.0 

7 

13.5 1. 

183346E+07 0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


FUEL . RLN file for First Sample Run: 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


FUEL Split: 3-1/ 1-0 
2 

4 . 055000E+03 486.0 1.185883E+07 

4 . 055000E+03 486.0 1.185883E+07 


2 

1 0 
8 2 


1 


17.97 

1.04154 0.0 

0.0 

0 


1.2785 

75.0 1.04154 

0.0 

1 


8.138 

1.04154 0.0 

0.0 

0 


1.2785 

-75.0 1.04154 

0.0 

1 


32.51 

1.04154 0.0 

0.0 

0 


1.2785 

75.0 1.04154 

0.0 

1 


8.65 

1.04154 0.0 

0.0 

0 


1.2785 

-75.0 1.04154 

0.0 

2 

1 

1 



1 


6.2 

0.52077 0.0 

0.0 

7 


4.5 1. 

183346E+07 0.0 

0.0 

2 

1 

2 



1 


6.2 

0.73648 0.0 

0.0 

7 


4.5 1. 

18334 6E+07 0.0 

0.0 

2 

3 




10 

0 



1 


17.97 

1.04154 0.0 

0.0 

0 


1.2785 

75.0 1.04154 

0.0 

1 


8.138 

1.04154 0.0 

0.0 

0 


1.2785 

-75.0 1.04154 

0.0 

1 


32.51 

1.04154 0.0 

0.0 

0 


1.2785 

75.0 1.04154 

0.0 

1 


8.65 

1.04154 0.0 

0.0 

0 


1.2785 

-75.0 1.04154 

0.0 

1 


6.2 

1.04154 0.0 

0.0 

7 


4.5 1. 

183346E+07 0.0 

0.0 


72.13 

72.13 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


II-2 



Input Unique to NYQUIST for First Sample Run 

CONST. RLN file for First Sample Run: 

0.1 6219.0 2.67 2.330000E-03 -315.0 

0.1 6219.0 2.67 2.330000E-03 -315.0 

0.1 6219.0 2.67 2.330000E-03 -315.0 

Input Unique to SSFREQ for First Sample Run 

IMODE.RLN file for First Sample Run: 


Sample Run 




2 

0.000000 

4.00000 

0.000697 

3.214000 

1.200000 

-315.0000 

450204.00 

450204.00 
0.01 

2.232000 

1716.000 
0.010000 

4000.000 

4000.000 

0.000000 

4.000000 

142500.0 

0.440000 

2.670000 

1965.000 

2 

0.000000 

4.00000 

0.000697 

3.214000 

1.200000 

-315.0000 

450204 . 00 

450204.00 
0.01 

2.232000 

1716.000 
0.010000 

4000.000 

4000.000 

0.000000 

4.000000 

142500.0 
0.440000 

2.670000 

1965.000 

2 

0.000000 

4.00000 

0.000697 

3.214000 

1.200000 

-315.0000 

450204.00 

450204.00 
0.01 

2.232000 

1716.000 
0.010000 

4000.000 

4000.000 

0.000000 

4.000000 

142500.0 
0.440000 

2.670000 

1965.000 


IHODE.FRQ file for First Sample Run: 
3 

110 

113 

116 

IMODE.TAU file for First Sample Run: 


0.0005 0.0007 0.0009 0.0011 0.0013 0.0015 0.0017 0.0019 
0.0023 0.0025 


.0021 
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II. 2 Files for Second Sample Run 
Common Input for Second Sample Run 
CLASS. ENG file for Second Sample Run: 


3.112000E+03 9.504000E+04 4.464000E+04 

CLASS. LOX file for Second Sample Run: 

Basic configuration 

4.055000E+03 2.264000E+03 1.185883E+07 7.213000E+01 

1 

1 1 

1 2.000000E+01 1.000000E+00 .000000E+00 .000000E+00 0.0 

1 2 . 000000E+00 1.000000E+00 .000000E+00 0.000000E-00 0.0 

1 4.000000E+00 1.000000E+00 .OOOOOOE+OO .000000E+00 0.0 

1 2. 000000E+00 1. OOOOOOE+OO .000000E+00 O.OOOOOOE-OO 0.0 

1 3 . 000000E+01 1. 000000E+00 .OOOOOOE+OO .000000E+00 0.0 

1 2. 000000E+00 1. OOOOOOE+OO .000000E+00 O.OOOOOOE-OO 0.0 

1 4. 000000E+00 1. 000000E+00 .OOOOOOE+OO .OOOOOOE+OO 0.0 

1 2. OOOOOOE+OO 1. OOOOOOE+OO .OOOOOOE+OO O.OOOOOOE-OO 0.0 

1 2.000000E+01 1. OOOOOOE+OO .OOOOOOE+OO .OOOOOOE+OO 0.0 

7 4 . 500000E+00 1. 183346E+07 0.0 0.0 0.0 

CLASS. FUL file for Second Sample Run: 

Basic Configuration 

4.055000E+03 8.48000E+02 1.185883E+07 7.213000E+01 

1 

1 1 

1 J ° 2 . 000000E+01 1. OOOOOOE+OO .OOOOOOE+OO .OOOOOOE+OO 0.0 

1 2. OOOOOOE+OO 1. OOOOOOE+OO .OOOOOOE+OO O.OOOOOOE-OO 0.0 

1 4. OOOOOOE+OO 1. OOOOOOE+OO .OOOOOOE+OO .OOOOOOE+OO 0.0 

1 2. OOOOOOE+OO 1. OOOOOOE+OO .OOOOOOE+OO O.OOOOOOE-OO 0.0 

1 3 . 000000E+01 1. OOOOOOE+OO .OOOOOOE+OO .OOOOOOE+OO 0.0 

1 2. OOOOOOE+OO 1. OOOOOOE+OO .OOOOOOE+OO O.OOOOOOE-OO 0.0 

1 4. OOOOOOE+OO 1. OOOOOOE+OO .OOOOOOE+OO .OOOOOOE+OO 0.0 

1 2. OOOOOOE+OO 1. OOOOOOE+OO .OOOOOOE+OO O.OOOOOOE-OO 0.0 

1 2 . OOOOOOE+Ol 1. OOOOOOE+OO .OOOOOOE+OO .OOOOOOE+OO 0.0 

7 4 . 500000E+00 1.183346E+07 0.0 0.0 0.0 

Input Unique to NYQUIST for Second Sample Run 

CLASS. CON file for Second Sample Run: 

1 . 000000E-03 6219.000000 2.670000 2.330000E-03 -315.000000 
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Input Unique to SSFREQ for Second Sample Run 
CLASS. IMO file for Second Sample Run: 


Basic configuration 
2 


0.000000 

4.00000 

0.000697 

3.214000 

1.200000 

- 315.0000 


95040.00 

95040.00 

0.01 

2.232000 

1716.000 
0.010000 


4000.000 

4000.000 

0.000000 

4.000000 

142500.0 
0.440000 


CLASS. FRQ file for Second Sample Run: 


2.670000 

1965.000 


3 

110 113 116 

CLASS. TAU file for Second Sample Run: 


0.0005 0.0007 0.0009 0.0011 0.0013 0.0015 
0.0017 0.0019 0.0021 0.0023 0.0025 



Appendix III 


Homework Problem: 

One configuration is given for all three codes, therefore the engine 
file and piping files need be developed only once. This is enoug 
for ADMIT, but NYQUIST and SSFREQ will require their unique files to 
be developed. The configuration consists of two LOX tanks, one FUEL 
tank and two engines. The two LOX tanks and the two engines are 
identical, one LOX line going to one engine. The fuel line has a 
split and feeds both engines. 

Assignment # l 

Develop the files for the engine (HW.ENG), the LOX line (HW.LOX), 
and the FUEL line (HW.FUL). Run both the LOX line and the FUEL 
line using ADMIT. Use a frequency range of 1 40 40. 

Assignment # 2 

Develop the file unique to NYQUIST and run the code using both 
FUEL and LOX lines. Use a frequency range of 1 40 300. 

Assignment # 3 

Develop the file unique to SSFREQ and run the code using files 
IMODE.FRQ and IMODE.TAU for the frequency and tau ranges. 

Description of Configuration 

LOX data (both tanks have the same values) - 
32 ■ diameter sphere 
mass flow rate of 3,000 lbm/sec 
density of 71.4 lbm/ft~3 
bulk modulus of 11,800,000 lbf/ft~2 
manifold volume of 15 ft~3 

manifold bulk modulus of 11,800,000 lbf/ft 2 
di am eter of piping from tank to engine is l ft 

FUEL data - 

two 32' hemispheres separated by 25' cylindrical section 

mass flow rate of 1,000 lbm/sec 

density of 72.1 lbm/ft~3 

bulk modulus of 11,800,000 lbf/ft~2 

manifold volume of 5 ft~3 

manifold bulk modulus of 11,800,000 lbf/ft~2 
diameter of piping from tank to split is 1.4 ft 
diameter of piping from split to engine is l ft 
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Engine data (both engines have the same values) - 

total mass flow of 3/500 lbm/sec 
chamber pressure of 660 lbf/in 2 
pressure drop across orifice of 312 • 5 lbf/in 2 

Combustion data - 

transport lag of 0.01 sec 
characteristic velocity of 6/200 ft/sec 
mixture ratio of 3 

characteristic time constant of .0023 sec 

change in velocity with mixture ratio of -300 ft/sec 

Chamber and Nozzle data (both engines have the same values) - 

length of combustion zone is 5 ft 
pressure is constant at 660 lbf/in 2 

temperature varies linearily from ejector (4,000 °R) to the 

nozzle (5,000 °R) . (Use 3 stations at x=0', 2.5', and 5.0 .) 
ratio of specific heats is 1.22 
chamber diameter of 3.2 ft 
throat diameter of 2.2 ft 
length of chamber of 5 ft 
invariant time lag of 0.0005 sec 
mixture ratio interaction index of 0.01 
gas constant of 1,716 (ft/sec) *2/°R 
maximum overpressure of 25,000 lbf/ft~2 
dc*/dr of -300 ft/sec 
dh/dr of 0.01 (ft/sec) ~2 

mass of liquid/unit chamber volume of 0.45 lbm/ft 3 
axial component of liquid velocity of 2000 ft/sec 
no damping 

The physical layout is shown in the figure on the following page. 




